Decay channels and charmonium mass-shifts 
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The discovery in the last few years of the X, Y and Z states of the extended charmonium family 
has highlighted the importance of the closeness of decay channels to an understanding of these 
mesons. We aid this debate by illustrating a simple calculational procedure for including the effect 
of open and nearby closed channels. 
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I. MODELLING DECAY CHANNELS 

The discovery of narrow states of hidden charm, the 
X, Y, Z- mesons [3, Q, has generated a whole literature 
discussing their nature, structure and relation to char- 
monium. The fact that a state, like the X(3872) sits 
between D*°W and D*+D~ thresholds Q, with a width 
of less than 1.2 MeV, has highlighted the potentially im- 
portant role that hadronic decay channels may have on 
the spectrum. Indeed, it is a feature of resonances with 
strong 5-wave thresholds that the states are drawn close 
to their strongly coupled thresholds [4[ as often discussed 
for the /o and ao close to KK threshold Q. Eichten, 
Lane and Quigg (ELQ) Q have calculated the effect of 
open channels for states with hidden charm in a scheme 
that many find unfamiliar. In this note we want to revisit 
an approach related to the Dyson summation for the in- 
verse meson propagator. This idea is not new and was 
considered for charmonium many years ago by Heikkila, 
et al. 7]. What is new here is the straightforward way 
in which we can estimate the effects of open and nearby 
closed channels. 

The inverse boson propagator, V(s), is shown in Fig. 1, 
where s is the square of the momentum carried by the 
propagator. With 11(a) the contribution of hadron loops, 
the complex mass function M(s) is related to this by 



V(s) = M 2 (s) - s = mf - s + n(s) 



71=1 



n„(s), (i) 



where uiq is the mass of the bare state and the sum is 
over all loops (Fig. 1). The propagator, V~ 1 (s), will then 
have a pole at (at least one) complex value of s = sr. 
This position specifies the mass and width of the phys- 
ical resonance. If we denote the threshold for the nth 
channel by s — s n , then clearly only those that are open 
for s ~ Re(s^) > s n contribute to the decay width of 
the physical hadron. However, in principle all hadronic 
channels contribute to its mass. Indeed, each of the in- 
finity of closed channels contributes to the real part of 
II(s) and for a given physical mass can be thought of as 
redefining the "bare" mass. Since we are interested only 
in mass-shifts, let us subtract Eq. (1) at some suitable 



point s = sq to be defined below, then 
M 2 (s) - M 2 (s ) = U(s) - n(a ) 



J2 P„W - n„( So )] 



(2) 



Since II ra (s) will be effectively constant for those virtual 
channels for which Re(sfl) << s n , their contribution 
will cancel out in Eq. (2). Consequently, the mass-shift 
is entirely given by the hadronic channels that are fully 
open or only just "virtual" . It is a reasonable expecta- 
tion that deeply bound states, like the J/ip, have masses 
defined by the charmonium potential. The mass of the 
J/4 1 then essentially defines the mass scale and fixes the 
charm quark mass at the relevant scale. It is thus natural 
to set so — M(J/ip) 2 . In line with expectation our re- 
sults change little if we use s — ^m 2 instead. For each 
state we take the value of A4(sq) to be that predicted by 
a charmonium potential, unperturbed by hadronic chan- 
nels. Of course, if the parameters in the charmonium 
potential are fixed with reference to physical states for 
which open charm channels may contribute, we have an 
issue of double counting. We believe that by fixing the 
charmonium parameters by only deeply bound states we 
avoid this problem. 

Since each II„ (s) is an analytic functions with a right- 
hand cut, we can write a Cauchy representation in sub- 
tracted form, so that 



An„(s,s ) = n„(s) - n„(s ) 



(s-s ) 



(s' - s) {s> - s ) 



Then 

E An «(vo) = M 2 (s)-Mr h 

n=l 



'! rmon ! a in 



= AM 2 (s). 



(4) 
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FIG. 1: The bare bound state propagator is dressed by 
hadronic loops. The dot signifies the dressed propagator and 
vertices. 
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The form of ImTI n for particle P coupling to each channel 
AB is taken to have a simple form, for s > s n : 

/2k\ 2L+1 

Imn„(s) = - gl \ -j=\ exp(-afc 2 ) , (5) 

where g n is the coupling of particle P to channel n (i.eto 
particles A and B), L is the orbital angular momentum 
between A and B, while k is the 3-momentum of A and 
B in the rest frame of P. So as usual 

Ak 2 /s = 1 - 2(m\+m 2 B )/s + {m 2 A - m 2 B f / s 2 . (6) 

The scale factor a is related to the radius of interaction, 
R, by a = R 2 /Q. This is in turn related to the size of 
the overlap between the cc and the AB states. A larger 
value of a produces a smaller mass shift. A value of 
a = 0.4 GcV~ 2 is favoured solely because it gives the 
most sensible results. This corresponds to R ~ 0.3 fm. 

For open channels, the coupling g n is simply related 
to the channel n decay-width through Eq. (3) with s ~ 
sr. For nearby closed channels we use the coupling to 
states with the same quantum numbers. As a guide to 
the size of the effects, the calculations presented here 
systematically include the channels DD, DD , D*D , 
and D S D S . 

II. COMPARISON WITH EXPERIMENT AND 
OTHER STUDIES 

We define the base from which the shift due to decay 
channels is to be computed by a non-relativistic poten- 
tial model for charmonium. From the many potential 
modellings we choose the classic work of Godfrey and 
Isgur [8[, more recently tabulated by Barnes, Godfrey 
and Swanson (BGS) [9]. This is because BGS not only 
provide a prediction for the eigenstates, but include cal- 
culations using the 3 Pq model for the partial widths. It 
is these that fix the couplings g n , which are the essen- 
tial input into Eq. (5) for computing the mass-shift from 
decay channels. Using these inputs we compute the cor- 
rection to the real and imaginary parts of AM 2 as shown 
in Fig. 2 for the example of the ifi'". From such plots we 
arrive at the mass-shifts given in Table I. These are pre- 
sented in two ways. The simplest is the shift in what we 
call the Breit-Wigner mass, for which we only need to 
compute An(s) at s = ReM 2 (s), Fig. 2. However, the 
physically relevant quantity is the shift in the position 
of the pole in the complex energy plane. This requires 
we evaluate AII(s) at s = sr = M. 2 (sr). For states 
with small widths, of course, the Breit-Wigner and pole 
masses differ little. However, for states with larger cou- 
plings, the difference is inevitably bigger. Indeed, some 
states get shifted below their threshold and their pole 
moves to the real axis. Others, however, are subject to 
significant changes. The largest effect is found for the 
X^, 1 (3 3 Pi) state, where a shift of Aitibw = —66 MeV is 
reduced to just Am po i e — —29 MeV. 
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FIG. 2: The real and imaginary parts of AM 2 (s) as functions 
of E = yfs for the ip"'(3 3 Si) propagator. The dashed line 
shows the curve — s, where it intersects the real part of 
AM 2 {s) labelled BW, defines the Breit-Wigner mass. The 
cusps in the real and imaginary parts occur at each of the 
thresholds. 



Correcting the bare masses delivered by the potential 
model of Godfrey and Isgur [8, 9] by our calculated decay 
channel induced mass-shifts brings better agreement with 
experiment as seen from Table I. For the the down- 
ward shift by between 36 and 41 MeV is reasonable. That 
for the 77" of 45 to 58 MeV is not quite enough to bring it 
in line with the measured mass, which is 100 MeV below 
the potential model prediction. 

States with common J PC quantum numbers have com- 
mon decay channels and so inevitably mix through these 
hadronic intermediate states. Two such states are the rj' c 
and rj". Since the ground state rj c is deeply bound, it 
mixes little with these. Explicit calculation gives a shift 
of -0.6 MeV for the rf c and even less for the r)". Con- 
sequently, for the states listed in Tables I and II, these 
inter-state mixings are small and can be neglected. 

We first compare our calculation with that of Heikkila, 
et al. Q of more than twenty years ago. These authors 
consider the spectrum of heavy quarkonia, in which the 
loop effects are built in from the start in the determina- 
tion of the parameters of the underlying non-relativistic 
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Name 


State 


Experimental Mass 


Potential Mass 


Thadrons 


ArriBW 


Ampole 




n 2S+1 Lj 


MeV 


MeV 


MeV 


MeV 


MeV 


Vc 


I 1 So 


2980±1 


2982 


- 


- 


- 


J/V 


l 3 Si 


3096.9 


3090 


- 


- 


- 




2 1 S 


3638±4 


3630 


- 


-10 


-10 


V.' 


2 3 Si 


3686.1 


3672 


- 


-9 


-9 


h c 


2 1 P 1 


3525.9 


3516 


- 


-2 


-2 


Xc 


2 3 P 


3414.8 


3424 


- 


-9 


-9 


Xci 


2 3 Pi 


3510.7 


3505 


- 


-16 


-16 




2 3 P 2 


3556.2 


3556 


- 


-6 


-6 


€ 


3 X S 


3943±6 


4043 


80° 


-45 


-58 




3 3 Si 


4039±1 


4072 


80±10 6 


-36 


-41 




3'Pi 


- 


3934 


87 a 


-5 


-12 




3 3 P 


- 


3852 


30 a 


-70 


-70 




3 3 Pi 


- 


3925 


168 a 


-66 


-29 




3 3 P 2 


- 


3972 


80° 


-55 


-48 




3 1 D 2 




3799 








r 


3 3 Di 


3771±2 


3785 


23±3 b 


-40 


-40 




3 3 D 2 




3800 










3 3 D 3 




3806 









TABLE I: Our Results. The calculated shifts for both the Breit-Wigner and pole masses are computed from a base defined by 
the non-relativistic model of Barnes, Godfrey and Swanson [sj listed above. Far from their masses, these predictions may be 
incorrect. Experimental data are from PDG Q|. Only experimental errors greater than 1 MeV are quoted. 
a BGS 0, b PDG @|. 



State 


Centroid 
MeV 


Spin Splitting 
MeV 


Bare Mass 
MeV 


Thadrons 

MeV 


ArriELQ 
MeV 


Our Mass 
MeV 


Our Am 
MeV 


2 1 So 
2 3 S'i 


3674 
3674 


-50.1 
16.7 


3623.9 
3690.7 




15.7 
-5.2 


3617.0 
3676.5 


-6.9 
-14.2 


&S 
3 3 S! 


4015 
4015 


-66 
22 


3949 
4037 


74 
49.8 


-3.1 
1.0 


3924.5 
4020.0 


-24.5 
-17.0 


3 1 Pi 
3 3 P 
3 3 Pi 
3 3 P 2 


3922 
3922 
3922 
3922 



-90 

-8 
25 


3922 
3832 
3914 
3947 


59.8 
61.5 
81 
28.6 


-5.4 
27.9 
6.7 
-9.6 


3892.0 
3818.8 
3868.9 
3939.4 


-30.0 
-13.2 
-45.1 
-7.6 


3 1 D 2 
3 3 P>i 
3 3 D 2 
3 3 D 3 


3815 
3815 
3815 
3815 



-40 


20 


3815 
3775 
3815 
3835 


1.7 a 
20. r 
0.045 
0.86 a 


4.2 
-39.9 
-2.7 
19.0 


3813.3 
3728.1 
3815.0 
3833.1 


-1.7 
-46.9 
0.0 
-1.9 



TABLE II: Comparison of the calculation and modelling by Eichten et al. in columns 2-5, with the results from our loop 
calculations from their same base bare masses with their channel couplings in columns 7, 8. The couplings to individual 
channels are taken from the partial decay widths computed by ELQ in their Table V of p]. The numbers with the superscript 
a are from Eichten et al. [lol ]. 



potential model. Meson loops then have a dramatic ef- 
fect on the "bare" states with shifts of hundreds of MeV 
in mass for the lightest states. In their calculation the 
infinity of virtual channels (or as many of these as they 
choose to include) all have an effect. In contrast, in our 



calculation by using subtracted dispersion relations for 
the meson loops, the effect of the many closed channels 
is absorbed into the subtraction constants. Moreover, 
because we expect deeply bound states like the J/ip to 
be negligibly affected by loop corrections and well ap- 
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proximated by charmonium potential calculations, the 
subtraction constants are accurately determined. The 
predicted widths by Heikkila et al. are within a factor 2 
of experiment for the ip" and ip'". 

We now compare our results with those obtained by 
Eichten, Lane and Quigg @. The first of two com- 
parisons is directly with their results given in Table II, 
columns 2-6. They ascribe part of the mass-shifts to spin 
splittings by suitably adjusting a s . With this the overall 
scale of their mass shifts assigned to decay channels is 
typically smaller — compare our results in the final col- 
umn of Table I with column 6 of Table II. For the t/j", 
we and they find a downward shift of 40 MeV. However, 
the most noticeable difference is that, in our modelling, 
loops shift the mass downwards, whereas Eichten et al. 
also have appreciable upward shifts. Though the decay 
rates for the 3 Pq and C 3 models are qualitatively similar, 
there are some important differences in mass-shifts. For 
example, for the 77", the couplings are computed from 
decay rates at different masses. For the BGS model, this 
is at a "bare" mass of 4043 MeV, when both DD* and 
D* D* channels f|[ are open. For the ELQ model, the 
mass is below D D* threshold. Consequently, the shift 
in pole position for the 77" is —58 MeV in Table I, and 
just —24 MeV in Table II. For the rj' c , the same chan- 



nels contribute, but both are virtual, and so we have 
assumed the couplings to these are as computed for the 
77". This results in a 10 MeV downward shift for the rj' c 
using the BGS couplings, whilst only a 2 MeV shift with 
the ELQ C 3 . This would reduce the spin-spin splitting 
by 7 or 8 MeV. A greater reduction is required. The 
experimental value is 48 MeV, whilst potential models 
typically predict between 60 and 80 MeV. 

A like for like comparison is to recompute the effect 
of meson loops in our calculational scheme but using the 
ELQ couplings and ELQ bare masses. The results are 
shown in Table II. Then columns 6 and 8 can be directly 
compared. We see a distinctly different pattern of mass- 
shifts with essentially only that for the ip" being similar. 
However, the calculation presented here is more straight- 
forward to reproduce and adjust to new information on 
partial decay rates to be measured in the future and so 
may serve as a simple guide to the size of decay channel 
effects. 
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